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We study the local flux of electrons and positrons from annihilating Dark Matter (DM), and in- 
vestigate how its spectrum depends on the choice of DM model and inhomogeneities in the DM 
distribution. Below a cutoff energy, the flux is expected to have a universal power-law form with an 
index n ~ —2. The cutoff energy and the behavior of the flux near the cutoff is model dependent. 
The dependence on the DM host halo profile may be significant at energies E < 100 GeV and 
leads to softening of the fiux, n < —2. There may be additional features at high energies due to 
the presence of local clumps of DM, especially for models in which the Sommerfeld effect boosts 
subhalo luminosities. In general, the flux from a nearby clump gives rise to a harder spectrum of 
electrons and positrons, with an index n > —2. Using the Via Lactea //simulation, we estimate the 
probability of such subhalo effects in a generic Sommerfeld-enhanced model to be at least 4%, and 
possibly as high as 15% if subhalos below the simulation's resolution limit are accounted for. We 
discuss the consequences of these results for the interpretation of the ATIC, PAMELA, HESS, and 
Fermi data, as well as for future experiments. 

PACS numbers: 95.35.+d, 96.50.S-, 98.70.Sa 
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I. INTRODUCTION 



the shape of DM host halo, and 



One of the corner stones of the standard big bang 
cosmology is the presence of dark matter (DM). Al- 
though DM comprises at least 80% of all matter in the 
Universe, surprisingly little is known about its nature. 
A recent rise of interest in the DM problem has been 
triggered by observations of anomalies in several high 
energy cosmic ray experiments. The PAMELA Collab- 
oration [l[ reported an increasing positron fraction in 
the energy range 10 - 100 GeV, while the ATIC ^, 
PPB-BETS 3, and HESS 0] experiments suggest a 
bump in the total e+e" flux between 100 GeV and 1 
TeV. Recent Fermi/LAT and HESS data show 
a smooth spectrum which is harder than the expected 
background at energies E > IQO GeV and has a break 
around 1 TeV but without significant bumps. 

The possibility that the anomalies are due to annihi- 
lating or decaying DM is very exciting, although stan- 
dard astrophysical sources, such as pulsars or nearby 

supernova remnants 3 B S, US Hil, Q Hi, Q, E H, 

provide a viable explanation as well. 

In order for DM annihilation to produce the neces- 
sary e+e^ flux, the cross section must be 100 - 1000 
times larger than the standard value {{crv) ^ 3 x 10~^^ 
cm^ s~^) inferred from the relic DM abundance. Boosts 
due to DM inhomogeneities are expected to be < 10 
[l3] in the solar neighborhood, and an additional mech- 
anism, such as the Sommerfeld enhancement, is neces- 
sary UM^- The Sommerfeld enhancement can be 
further amplified in subhalos due to their lower veloc- 
ity dispersions [2]|. DM models, proposed to explain 
the ATIC and PAMELA anomalies, include the light- 
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25 



26 



est supersymmetric particles 

Kaluza-Klein particles |28l. l29l. 30 1, and various phe 



22, 23, 



27| 



33|,lM,l35|,l36|,l37| (for 



nomenological scenarios |31|. I3S 
a recent review see, e.g., [38|, |39|) 

In the presence of a large number of competing pos- 
sibilities, one of the most important questions is what 
properties of the e+e" flux are model dependent and 
what properties are universal. In this work we study 
the dependence of the e+e^ flux on 

• DM models. 



• the presence of subhalos. 

The dependence on DM models and on the host halo 
density proflle was recently investigated in 



31, 40 



41 



421 . One may notice that the fiux is nearly independent 
of the density profile and that the choice of DM model 
seems to affect the spectral shape only at high energies 
near the break. 

We show that for any DM model there exists an 
energy E^, < A/dm such that the e+e^ flux has 
a universal behavior Fe± (E) ^ E^"^ for energies 
1 GeV <^ E <^ E^. In general, for E <^ E^, the index 
depends on the energy loss function, E = —b{E), for 
electrons propagating in the interstellar medium (ISM) 



HE) 



E<^E^ 



(1) 



where h{E) E'^ for > 1 GeV. The behavior in Eq. 
(HI) does not depend on the DM model, but the value of 
E^ is model dependent. Eq. ^ may be also modified 
at energies E '-^ IQ GeV due to a dependence on the 
DM density profile. 

This paper is organized as follows: In Sec. |TT1 we 
review the propagation of electrons and positrons in 
the ISM and derive Eq. In Sec. Iml we study 

the dependence of the spectral index on DM models 
and on the spatial distribution of the DM. For our 
purposes, the most important difference between DM 
models is the number of steps in the annihilation-decay 
process, with more decay steps leading to a softer e^e~ 
injection spectrum. However, propagation effects then 
ensure that the final spectrum exhibits the universal 
behavior of Eq. ([T]). Likewise, variations in the DM 
spatial distribution also do not change this behavior. 
As we show, the spectral slope is not affected by lin- 
ear gradients in the source function, with the leading 
order correction arising from second derivatives. The 
Earth is sufficiently far from the Galactic center such 
that the corresponding effects are small, although the 
corrections may become important at £' ^ 10 GeV due 
to an increase in the characteristic diffusion distance. 

A local clump of DM, on the other hand, may re- 
sult in a significant variation in the spectral index. 
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In Sec. IIV( we investigate this possibility using Via 
Lactea II (hereafter VL2), one of the highest resolution 
numerical simulations of the DM structure of a Milky- 
Way-scale halo 17|, |43|, |44| . We find that on average 
there is one M > IO^Mq clump within 3 kpc from 
Earth. Without Sommerfeld enhancement, the likeli- 
hood of boosting the local flux by a significant factor 
(>10) due to one such nearby subhalo is very low, in 



In the examples, we will use Dq 

5 



0.6 



M ISa, l54| 



20 pc^ kyr ^ and 
We will also assume that the 



agreement with previous findings [17|, |45| . Allowing for 
Sommerfeld enhancement, however, the average lumi- 
nosity from nearby clumps can match (or even exceed) 
the smooth host halo contribution. We show that, in 
general, the electron flux from a nearby DM clump is 
harder than from the host halo (as also observed in cording to the diffusion- loss equation (49|, |50 1 



diffusion zone in the vertical direction (away from the 
Galactic plane) is sufficiently large, L > 4 kpc, such 
that the escape of electrons from the Galaxy can be 
neglected. Models with a small diffusion zone, L < 1 
kpc, tend to have exponential suppression of the flux 
from DM at energies < 100 GeV which would 
contradict the PAMELA data. Our main qualitative 
conclusions will not depend on the choice of parameters 
bo, Dq, and 6. As a useful simplification, we will also 
assume that these parameters do not depend on x. 
The electron energy spectrum / — evolves ac- 



dS-4 



45 
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48|). In the DM annihilation picture, the 
high energy (- 600 GeV) flux observed by ATIC, PPB- 
BETS, and HESS is consistent with being produced by 
a single DM clump, but in order to explain PAMELA 
the flux must be dominated by the smooth host halo 
below energies of ~ 100 GeV (another possibility to 
explain PAMELA is via a second, more distant, large 
clump of DM as pointed out in [45] ) . Recent Fermi and 
HESS data are significantly smoother than the ATIC 
data and are fitted better with the flux from the host 
halo without contributions from large subhalos. We 
conclude in Sec. |V] by summarizing which properties 
of the e^e^ flux are model dependent and which are 
universal. 

II. PROPAGATION OF ELECTRONS 

The propagation of electrons in the ISM is described 
by diffusion in the Galactic magnetic field and energy 
losses 



49 



50| . At energies E > 10 GeV, the main losses 
are from synchrotron radiation and inverse Compton 
scattering, which can be estimated as 

E = -baE\ (2) 

where 6o = 1-6 x 10^^^ GeV^^s^^ for the local densi- 
ties of starlight, infrared, and CMB photons [5l| and 
a magnetic field of 3^G [50|. For relativistic electrons, 
the diffusion coefficient can be expressed as 

DiE)^DJ-^)\ (3) 



a/ 

dt 



d_ 

'dE 



{h{E)f)+D{E)V^f + Q{^,E,t). (4) 



The general solution of this equation is found in 
[49I . [sq . If the source is constant in time, then 
df /dt = and the equation simplifies to the usual dif- 
fusion equation. If we introduce a new variable \{E) 
such that 



then 



d\ _ D{E) 
dE ^ ^ b{E) ' 



dX ' D{E) 

The Green's function for the diffusion equation is 



G(x,A) 



1 



-e 1^ . 



(5) 



(6) 



(7) 



(47rA)3/2 

and so the general solution to Eq. ^ is given by 

/(x,^) = ^Jd^xoJ^ dEoG{x-xo,HE,Eo)) 
•Q(xo,So), (8) 



where 



X{E,Eo) 



IE 



^° D{E')dE' 
KE') ' 



(9) 



Our approximation is valid if the characteristic propa- 
gation distance, = 6X{E, £'o)|£;p^oo, is much smaller 
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than the size of the vertical diffusion zone L. For our 
choice of parameters, this gives 



S > 60 



,kpc 



GeV 



(10) 



If L ~ 4 kpc, then i? > 30 GeV. For smaller energies 
there is an exponential suppression due to escape of 
electrons from the Galaxy. 

For a homogeneous monochromatic source 
Q(xo,i?o) — Qo^iE — Mdm), the solution to Eq. 



fiE) = 



(11) 



which has the behavior ([T]) announced in the Introduc- 
tion. For energies S > 1 GeV, b{E) - E'^ (cf. Eq. ^) 
and therefore f{E) ^ E^'^ . One of the main purposes 
of the paper is to elucidate how the index n = — 2 
changes for different DM models and different DM dis- 
tributions. 



III. DARK MATTER IN A HOST HALO 

The source function of electrons and positrons com- 
ing from annihilating DM is 



1 Pdm 



dE 



(12) 



where m-^Ii number density, {av) is 

the thermally averaged annihilation cross-section, and 
dNann/dE IS the differential energy spectrum of elec- 
trons and positrons produced in a single annihilation 
event. Here we assume that the DM particle is its own 
antiparticle, e.g., a Majorana fermion, otherwise there 
is an additional factor of 1/2. 

As we discuss in more detail below, the cross sec- 
tion may exhibit a velocity dependent boost factor 
(av) — {av)o S{v), where (<7v)o is the cross section 
at freeze-out and lim S ^ 1. Generically, this enhance- 

V — *c 

ment scales as ^ 1/w, although resonances exist for 
certain parameter choices, in which case S ^ l/v"^ [l^. 
Note that in cither case the boost saturates at small 
velocities. 



DM halos in general are not isothermal and have 
lower velocity dispersions, and hence higher boost fac- 
tors, in their centers 



57 



58|. Given a density pro- 
file, and assuming equilibrium and spherical symmetry, 
it is straightforward to solve the Jeans' equation for the 
corresponding velocity dispersion profile crv(r). For a 
Navarro, Frenk, and White (NFW) density profile, the 
resulting velocity dispersions peak at about the scale 
radius and decrease towards the center as crv(r) ^ r^/^ 
[sot . This is a limiting case [i^l, and all other pro- 
files (e.g. the Einasto profile) have shallower velocity 
dispersion profiles. 

The source function (|12p can be split into a product 
of x-dependent and iJ-dependent functions 



where 



Q{x,E)^ kL{k)Q{E), 



9M^ 



(13) 



(14) 



QiE)^^{E), 



The x-dependent part is the luminosity 

L(x) = S{v{^)), (15) 

and the iJ-dependent part is the injection spectrum 

dE -'^ 

where we choose the normalization such that 
Q{E)dE is the total average number of electrons 
and positrons produced in an annihilation. 

A. Sommerfeld enhancement 

In the following we will assume the standard DM 
density pdm — 0.3 GeV cm~'^ and a freeze-out cross 
section {(jv)q = 3 x 10~^^ cm^ s^^ . This cross section 
is too small to give a sufficient annihilation rate for 
the PAMELA and ATIC anomalies, and hence a boost 
factor of 100 - lOOO is necessary. Recent Fermi data 
are consistent with this conclusion. If we assume that 
the PAMELA anomaly is due to dark matter, then the 
absence of steplike features in Fermi data below the 
cutoff requires Mdm ^ 1 TeV. Consequently the same 
boost factor 100 - 1000 is required to explain PAMELA. 
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An elegant way to increase the annihilation cross sec- 
tion without affecting the relic abundance, is to assume 
the existence of a new light force carrier (p in the dark 
sector, 7710 ~ C'(GeV). This additional force results in 
a Sommerfeld enhancement [l8| , which is unimportant 
at freeze-out when particles are close to relativistic, but 
can significantly boost annihilation rates today. There 
may be further enhancement, if the annihilation pro- 
ceeds through the formation of WIMPonium [6l|, a 
meta-stable bound state of two DM particles. This 
process is analogous to the annihilation of electrons 
and positrons through positronium, which is dominant 
at low velocities. For some parameters in the DM sec- 
tor, it is possible to radiatively create WIMPonium by 
emitting a particle, and in this case the WIMPonium 
annihilation channel dominates over immediate annihi- 
lation. 

We now demonstrate that in every model with pa- 
rameters that do not allow for WIMPonium creation, 
there exists an upper bound. 



(17) 



on the nonresonant boost factor, where v is the relative 
velocity of the particles. 

Generic (nonresonant) Sommerfeld enhancement is 

H 



Q( \ '^'^ 

V 



(18) 



This enhancement factor saturates when the deBroglie 
wavelength of the DM particles becomes equal to the 
force range ~ thus 



S < 



(19) 



WIMPonium cannot be created if the corresponding 
binding energy is smaller than the mass of the mediator 
1 2 

-a A/dm < niA,. 
4 ^ V 

Collecting these pieces together, we find 



(20) 



and thus 



na ^ TraA^DM ^ 47r 
V ^ ^ a ' 



(21) 



which gives the bound (I17p . In any particular model 
the precise value of the numeric coefficient in p7p may 
be different, but the parametric form should be the 
same (unless there are some additional reasons that 
prevent WIMPonium creation). 

At the location of the Sun (8 kpc), the velocity dis- 
persion of the DM particles is CTv ~ 200 km/s 2l|. 
Taking this as a proxy for the relative velocities of an- 
nihilating DM particles, we find 



S'(200 km/s) < 250. 



(23) 



Thus, the local Sommerfeld enhancement, 5(200 
km/s), cannot exceed ~ 250 without taking into ac- 
count WIMPonium creation or resonance effects. If we 
allow for WIMPonium creation (m^ < I/Ao^AIum), 
then the annihilation cross section increases by an ad- 
ditional factor of 7 (3) for fermionic (bosonic) DM par- 
ticles 6l|, i.e., for fermionic DM particles the enhance- 



ment factor is 



Siv) w 20- 



(24) 



In order to parameterize the Sommerfeld enhance- 
ment as a function of velocity, we need to specify 
the normalization and the saturation velocity Wmin- 
Then, for small velocities (w <C a), 



Siv) = So 



V + v„ 



(25) 



We take vq to be the local velocity dispersion, vq — 
200 km/s. In the examples we discuss below, we set 
v-min = 3 km/s and 5*0 = 600. These parameters are 
rather generic and can be obtained, for instance, in 
a fermionic DM model with a vector boson mediator 
with a — 1/50 and m^/M-oM = 10~^. The small mass 
ratio is necessary to have a small saturation velocity, 
^^^^ ^ ■ For small Wmin the role of small subhalos 

c Mum 

is greater since their velocity dispersion is smaller. For 
larger Wmin the luminosity of small subhalos becomes 
smaller and the probability to have an observable effect 
is reduced. 

For A/dm ^ 1 TeV we have ^ 10 MeV > 2m^± , 
i.e. this model is kinematically viable. In building a 
realistic model, one has to check various constraints 



(22) from diffuse gamma rays and neutrino fiuxes 4l|, |42 
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62, |63|, |6J] . However, instead of trying to find a specific 
model that satisfies all current observational data, we 
now discuss some general properties of the e+e^ flux 
from annihilating DM. 

B. DM model dependence 

For a homogeneous monochromatic source, the elec- 
tron spectrum was derived in Eq. (|f f p . In general, the 
spectrum can be approximated as f ^ E" with an en- 
ergy dependent index n{E). Since at low energies we 
expect f{E) ^ l/b{E), the variation of the index will 
be defined as 



_ d\ogb[E)f[E) 
dlogi^ ■ 

For a homogeneous source, the electron density is 



(26) 



I{E) = 



1 



Q{Eo)dE„. 



(27) 



biE) JE 

If the injection spectrum has the form of a power law 

Q{E) r^E" a> -1, (28) 

then the integral is saturated at Eq ^ M-qm and the 
variation of n is small for E <^ A/dm • 

We define the break energy E^, by the condition 



Sn{E^) = -1. 



(29) 



The flux Fe± ^ E ^ near E^, which is approximately 
the same scaling as the backgrounds 0, EHl ■ Thus the 

F ± 



ratio 



is maximal near E^ 



-fbackgr 

As an example, we consider DM annihilation fol- 
lowed by a chain of (fc — 1) two-body decays 

2x ^ 201 ^ 402 ^ ... ^ 2''-^e++2''-'^e- (30) 

We assume that rrii+i <C m^, i.e. the decay products of 
are relativistic in the rest frame of 4>i. The spectrum 
of 0i particles is [42 1 



UE) - 2SiM^~E) 

^ (log Mom ^ 



QriE) = 



i~2 



E 



i > 1, 



(31) 



where we include the multiplicity 2* in the definition of 
the source functions Qi{E). 
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FIG. 1: The e e~ flux depending on the number of steps 
in the toy model DM annihilation (see text). In order to 
keep the maximum values of E^F at _E* = 600 GeV, we 
choose the DM mass A/dm ~ e 2 E^- Apart from overall 
normalization, the spectra look similar to each other for 
E <^ Et. Model dependence shows up for energies E > E*. 



The DM particle mass is not directly observed, 
rather the experiments measure the peak of the bump, 
i.e. i?*, which is a function of Mdm. Motivated by the 
ATIC data we fix = 600 GeV and find the necessary 
Mdm in order to satisfy condition 



M] 



(32) 



= 1, 2, 3.5, 6, 10 ... 
k = 1, 2, 3, 4, 5 ... . 

For large k, the ratio grows exponentially, 
e . The corresponding fluxes for k < 5 are shown 
in Fig. [TJ At energies E <C 600 GeV the spectra look 
very similar to each other, indicating model indepen- 
dence. At higher energies, on the other hand, the spec- 
tral shape is sensitive to the value of k, i.e. it depends 
on the DM model. 



C. Profile dependence 

In this section we calculate the variation of the index 
due to a coordinate dependence of the luminosity. In 
order to find the spectrum of electrons coming from 
DM annihilating in the host halo, we take the source 
function (fT3|) and substitute it in Eq. ([8]). At the 
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observer's position x = 0, we get 



b{E)f{E) 



dEo kQ{Eo) 



(fxQ 



(47rA)3/: 



re-w i(xo) (33) 



In the previous section we considered the variation of 
the index due to the lower hmit of the Eq integration. 
Now the integral over xq is a function of X{E, Eq) and 
this introduces an additional E dependence. 
The resulting variation is 



Sn 



d\ogX (fxi)-3(L) 



(34) 



where 4^ 



dlogE (L) 
diogA _ „ ^^(j ^YiQ average of a function 



llogE 

f{x) is defined as 
(/> 



d^x 



(47rA)3/2 



e " f{x). 



(35) 



In particular, (x^) = 6A. Thus, for a uniform source, 
5n = 0. For a linear gradient 5L = bix'^ , (SL) = 
(^(5iy) = 0. Consequently, the first correction to n 
follows from the second derivative of L{x). 
For slowly varying i(x). 



(36) 



At the position of the Earth, the hmrinosity can be 
approximated as 



(37) 



Consequently 



V'Lir) = ^OyOr (pImS) - ^^^^L{r). (38) 



Substituting this expression in Eq. (|36|) . we get 

6n ~ —(1 — 6)a(a — 1)^, 



(39) 



where rg « 8 kpc is the distance from the center of 
the Galaxy to Earth. Since A <C rg, this variation is 
generally small (Fig. [2]). 

The situation is different when the density inhomo- 
geneity is due to a discrete subhalo. For a DM clump 
with size Id <^ \/A we can neglect the first term in the 
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FIG. 2: The dependence of e'^e flux on the DM halo pro- 
file. We use NFW [gi], Einasto {a = 0.17) [gi, and Isother- 
mal [gsI . [ggI profiles. The difference from the homogeneous 
distribution is proportional to aj^jg (i5)/ro, where ro = 8 
kpc is the distance from the center of the Galaxy. For small 
energies the characteristic diffusion distance increases and 
the corrections are more significant. Similar results were 
obtained in ^ for Mom = 300 GeV. 



numerator of Eq. ([34]). The correction to the index then 

becomes positive and 0(1), Sn « (1 - ^)- 

This leads to a harder electron spectrum, and the in- 
dex is no longer independent of the DM model and the 
properties of ISM. We consider some examples in the 
next section. 



IV. CLUMPS OF DARK MATTER 

In this section, we describe the properties of subhalos 
using the results of the Via Lactea II (VL2) simulation 
and study their influence on the electron flux spectrum. 

A. Via Lactea subhalos 

With a particle mass of 4000 Mq, VL2 can resolve 
subhalos down to ~ lO^M©. For our purposes, the 
important subhalo properties are their size, their num- 
ber density, and their luminosity. Since the annihila- 
tion rate scales as p^S, it is strongly peaked toward 
the centers of subhalos. For an NFW profile, ^ 90% of 
the total annihilation luminosity originates from within 
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FIG. 3; Some properties of DM subhalos using results of 
Via Lactea II simulation [l^. The upper plot shows the 
expected number of subhalos within a distance r from the 
Earth. In the lower plot we show the ratio of luminosities of 
all subhalos within r to the total luminosity of the host halo 
in the solid sphere of radius r around the Earth. We plot 
the ratio of luminosities without Sommerfeld enhancement, 
with generic S ^ 1/v enhancement, and with resonant S ~ 
enhancement. These curves only account for subhalos 
resolved in the Via Lactea //simulation, Msub ^ 10^ Mq, 
and should be viewed as lower limits. 



the scale radius, which is much smaller than the typi- 
cal propagation distance of ^ 1 kpc. Consequently, we 
will treat the subhalos as point sources. 

We determined the properties of the local subhalo 
population from the simulation by randomly placing 
100 sample spheres of radius 20 kpc, each at a dis- 
tance of 8 kpc from the host halo center. For each 
sample sphere, we identified all subhalos whose centers 
fall within the sphere. The top panel of Fig. [3] shows 
the mean, root mean square, minimum, and maximum 
number of subhalos over the ensemble, as a function of 
radius within the sample sphere. To a good approxi- 
mation the mean number of subhalos grows as , and 



hence we will assume that the subhalo number density 
is constant. On average there is one subhalo within 
r « 3 kpc. Note that these results refer solely to the 
portion of the substructure hierarchy that is resolved 
in the VL2 simulation, i.e. Aigub > 10^ Mq. The sub- 
halo mass function has been measured in simulations 
to be a simple power law dn/dM ~ M^" over 4-5 
decades of subhalo mass 
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71j . with a logarithmic 



slope a Ki 1.9 independent of distance from the host 
halo center. This implies that the local density of sub- 
halos with mass larger than M may be extrapolated to 
M < W^Mq as 



nsub(M) « 9 X 10"^ kpc 



-3f M 



V 1O5M0 



-0.9 



(40) 



In the bottom panel of Fig. [3] we plot the mean 
and the root mean square of isub/-^host) the ratio of 
the luminosity from all VL2-resolved subhalos to the 
smooth host halo luminosity, as a function of the en- 
closed radius within the sample spheres. We assigned 
a luminosity to each subhalo by assuming NFW den- 
sity and velocity dispersion profiles that are matched 
to the values of Vmax and rvmax as measured in the 
VL2 simulation. Vinax is the peak of the circular ve- 
locity curve v1{r) = GM{< r)/r, and rvmax is the 
radius at which this peak occurs, and both quantities 
are robustly determined in the numerical simulation. 
We have checked that our results do not change quali- 
tatively if an Einasto profile is assumed instead. 

We compare three different models: one with S — 1 
(no Sommerfeld enhancement), one with S' ~ 1/u, and 
one with S ^ Xjv^. Without Sommerfeld enhance- 
ment, the average contribution to the total luminosity 
from M > lO^AfQ subhalos is negligible, in agreement 
with previous findings In the S ^ Xjv model, 

however, this subhalo contribution is already as large 
as that from the smooth host halo, and in the S ^ l/v^ 
case it dominates by close to 2 orders of magnitude. We 
caution that these are ensemble averaged ratios: even 
at a distance of 3 kpc only 70% of the sample spheres 
contain one or more M > IO^Mq subhalo, and at 2 
(1) kpc this fraction drops to 24% (3%). The situa- 
tion for models with Sommerfeld enhancement can be 
summarized as follows: 
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i) when one or more M > 10^ AIq subhalos happen 
to he within the electron difFusion region, their 
combined flux dominates that of the smooth host 
halo; 

ii) the probability of finding one or more such sub- 
halos remains below 50% out to 2.6 kpc, and a 
mean occupancy of one such subhalo is reached 
only at 3 kpc; 

iii) the ensemble averaged expectation value of 
isub/-^host is unity for 1/v models (Af > IO^Mq 
subhalos are about as important as the host halo) 
and - 100 for l/w^ models (M > IO^A/q subha- 
los are much more important than the host halo) . 

These results suggest that if a model relies on the Som- 
merfeld effect to explain the high energy cosmic ray 
anomalies, then there is a non-negligible probability 
of an individual DM subhalo affecting the high energy 
electron flux. 

We stress once more that these results refer merely 
to the subhalos resolved in VL2. In the cold dark 
matter picture of structure formation, the clumping 
of DM should continue far beyond the artiflcial resolu- 
tion limit of current state-of-the-art numerical simula- 
tions like VL2. The total luminosity from all subhalos 
below the simulation's resolution limit may very well 
dominate the local pair production rate from DM an- 
nihilation, if the dense cores of such low mass subhalos 
are able to withstand the disruptive forces from tidal 
interactions with the DM host halo and with stars and 
molecular clouds in the Galactic disk '72]. However, 
the diffusive nature of the electron propagation would 
average out the contribution from subhalos with an 
0(1) probability of lying within the diffusion distance 
of '^l kpc, rendering their contribution indistinguish- 
able from the homogeneous host halo source. Let us 
define the characteristic mass scale Mi such that 

-Tr4m / dnsnh = 1. (41) 

«J J Ml 

where x^m ^ 1 kpc in our case. For the local sub- 
halo abundance of Eq. (gO]), Ah = 2,6OOM0. Any 
additional features in the electron spectrum from an- 
nihilating DM must arise from a small number of large 
subhalos, M > Mi, within the diffusion region. 




s 



FIG. 4: The probability to observe an order one feature 
at high energies from a clump of dark matter versus the 
index s of the luminosity function [Eq. (|42|l ]. The break at 
s — 2 corresponds to equipartition of DM luminosity. For 
s < 2 the luminosity is saturated by large clumps, for s > 2 
the luminosity is saturated by small clumps. The overall 
normalization is model dependent. 

B. The flux from a local clump 

In this section we study an example of the flux from 
a single DM subhalo in an S" ~ l/v model and compare 
it with the flux from the host halo. 

At flrst, let us estimate the probability to see a fea- 
ture from a single DM clump relative to the flux from 
the host halo and the averaged flux from minihalos (de- 
flned to have M < Mi ) . Effectively, the answer to this 
question depends on only one parameter, the logarith- 
mic scaling s of the subhalo number density with re- 
spect to the luminosity 



Qualitatively, if s > 2 then the integrated luminosity is 
saturated by small clumps and the probability to have 
a signiflcant flux from a large clump is negligible. If 
s < 2, then the luminosity is saturated by the large 
clumps and the probability to see an additional feature 
is significant. 

Let us estimate the probability of an 0(1) feature 
for the equipartition distribution, s — 2. In this case 
nsuh{L) ~ L~^, i.e. subhalos with 10 times larger lu- 
minosity have 10 times smaller density. Denote by Li 
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a subhalo is 




FIG. 5: The luminosity function from VL2 subhalos, as- 
suming NFW density and velocity dispersion profiles for 
the same three models as in Fig. |3] Without Sommerfeld 
enhancement (circles) the luminosity function is equipar- 
tioned (s = 2) down to the simulation's completeness limit, 
L ~ IO^Mq pc~^. For comparison we also show the distri- 
bution when the sample is restricted to subhalos with more 
than 250 particles (M > IQ^Mq) [small circles]. In the 
S ^ 1/v Sommerfeld case (squares), the luminosity func- 
tion is steeper (s = 2.8) above the saturation luminosity of 
Lsat ~ 2 X 10^ Mq pc~"^. Below saturation the distribution 
is expected to flatten to s = 2 (dashed lines), and indeed 
this behavior is clearly seen for a model with a higher satu- 
ration luminosity of Lsat ~ 3 x W^^A'Iq pc~'^ corresponding 
to m^/A/DM = 5 X 10~^ (small squares). For S ^ (tri- 
angles) the luminosity function above Lsat = 10^^ Mq pc~^ 
is even steeper (s = 4.9). 



the luminosity corresponding to a subhalo of mass Mi . 
Suppose that there are p decades of subhalos below 
Ml ~ 10"^ M0, which in the equipartition case each con- 
tribute an equal amount to the total luminosity. For a 
minimum DM clump mass of lO^^^ - IO^^'Mq |73l.l74| 
we can roughly estimate p — 7 — 15. 

The total luminosity of subhalos with L < Li is then 
-^minihaios = pLi- For the purposc of estimation, sup- 
pose that the luminosity of the host halo is smaller than 



iminihaios; then we should expect to observe an 0(1) 
feature if there is a subhalo with luminosity Lp — pLi 
within the diffusion distance. The probability of such 



^ ^ nsuhjLp) _ 1 

^ 71sub(-^l) P 



0.1, 



(43) 



and is independent of the value of Li and of the details 
of the propagation model. 

For general s, the luminosity of a subhalo that is able 
to produce a significant feature is 



Lp{s) 



;ub 



s-2 



Li. 



(44) 



The luminosity Li is defined as ^a^^jg duguh = 1 
and ccdiff 1 kpc in our case. The corresponding den- 
sity is 



and the probability becomes 



7r(s) 



nsuhjLp) _ / Lp{s) 
nsuh{Li) \ Li 



(45) 



(46) 



The normalization of this probability is model depen- 
dent. In particular, it depends on the ratio of the host 
halo luminosity to that in the subhalos. The shape of 
7r(s), however, is universal. It has a break at s = 2 
with a flat probability for luminosities saturated by 
large clumps, s < 2, and a rapidly decaying probabil- 
ity for luminosity saturated by small clumps, s > 2. 
The function 7r(s) for Li/L^i^ ^ 10^° is presented in 
Fig. HI 

In Fig. [5] we show an empirical determination of the 
subhalo luminosity function dNsnh/dL from the VL2 
subhalos for the same three models we considered in 
Fig. [21 As before, the subhalo luminosities are esti- 
mated using an NFW density and velocity dispersion 
profile, and we use the S ^ 1/v model normalized to 
S'o = 5'(200 km/s) = 600. We include all subhalos 
within the host's virial radius r2oo — 402 kpc. In the 
case without Sommerfeld enhancment we find s = 2 
(equipartition) down to the luminosity completeness 
limit of L ~ lO''^ Mg pc~'^. The turnover at smaller lu- 
minosities is due to the fact that we do not resolve 
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FIG. 6: The flux of electrons from a typical DM subhalo 
at distances r = 0.5, 1, 2 kpc from Earth versus the flux 
from the DM annihilation in the host halo. The flux from 
a subhalo is significant for r < 1 kpc and the corresponding 
spectrum is harder than the spectrum from the host halo. 
We use S ^ Ijv Sommerfeld enhancement of annihilation, 
Mdm = 2.1 TeV, m^/MoM = 10"^ a = 1/50, and 5'o = 
600. The subhalo luminosity is given in Eq. (|47|l . 



halos with M < IO^Mq and correspondingly low lu- 
minosities. When only subhalos with more than 250 
particles (M > IO^Mq) are included, then the distri- 
bution departs from a power law at roughly 10 times 
higher luminosities. 

The Sommerfeld-enhanced luminosity functions are 
steeper, with s = 2.8 for S ~ 1/v, and s — 4.9 for 
the S ~ 1/w^ model. In the smallest mass subhalos re- 
solved in VL2, the internal velocity dispersions become 
comparable to the saturation velocity. Even smaller 
subhalos will be fully saturated, and the slope of their 
luminosity function should be identical to the non- 
Sommerfeld-enhanced case. We thus expect a break 
in the power law right around the saturation luminos- 
ity Lsat- We estimate Lgat by determining the mean 
luminosity of all subhalos with T^nax/c within 10% of 
Wmin/c = to^/A/dm = 10"^, and find isat = 2 x 10^M|, 
pc-3 for S ^ 1/v and L^at = 10"M2 pc'^ for 
S" ~ 1/ti^ . Since these values are very close to the VL2 
completeness limit, it is difficult to clearly distinguish 
this break. For a model with a higher force carrier 
to DM particle mass ratio of m^/MoM = 5 x 10~^, 
this break is indeed apparent at the higher saturation 
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FIG. 7: Electron and positron flux from annihilating DM 
in the host halo plus a 1 kpc clump shown in Fig. |6] The 
hard spectrum of the flux from the subhalo is consistent 
with the ATIC bump around 600 GeV, but, in order to fit 
the PAMELA data, we need the softer fiux from the host 
halo at energies E < 100 GeV. We use the primary electron 
background ~ E~^'^ and secondary electron and positron 
backgrounds ~ E~^'^. The parameters of the DM model 
are the same as in Fig. |6l 



luminosity of Lgat = 2 x 10^M|, pc~'^. These results 
indicate that subhalos below the VL2 resolution limit 
are unlikely to either strongly enhance or suppress the 
likelihood of having a significant spectral feature due 
to an individual subhalo. 

As a concrete example, we consider the electron 
flux from DM annihilating in a clump at a distance 
a^sub = 1 kpc. We take the S ^ 1/v model with the 
normalization 5*0 = S'(200km/s) = 600 and choose the 
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FIG. 8: The same as in Fig. [7] but without the contribu- 
tion from the subhalo and with a larger local Sommerfeld 
enhancement in the host halo So — 1000 corresponding to 
a = 1/30. 



following luminosity of the DM clump 



(47) 



This luminosity is equal to the total luminosity in the 
host halo within approximately 2 kpc from Earth. Lgnh 
is about 4 times smaller than the luminosity resolved 
in VL2. Consequently, in the equipartition case s — 
2, the probability to find a clump of DM with L > 
Lsub is around 4 • 4% = 16%. The source function for 
the electrons and positrons from DM annihilating in a 
subhalo at position Xgub is 

Q,ub(x, E) = KL,„b<^'(x - x,ub)Q(^). (48) 

The flux from a DM subhalo can be obtained by sub- 



stituting the source Qsub into the general solution 
c 1 



sub 



47r b{E) 



kL 



sub 



dEo 



1 



(47rA)3/2 



QiEo). (49) 



In our example, we use the DM toy model from sec- 
tion nil Bl with the annihilation chain 2% — > 2(/)i 
A(j)2 — * 4e++4e^. The corresponding source function 
Q3{E) was presented in Eq. (PT|l . This model can be 
considered as a toy model for a more realistic decay 
through pions or muons instead of (j)2 ■ If there are four 
pions (muons) after the decay of 0i , then there will be 
only 2 electrons and 2 positrons in the final state and 
the boost factor in this model should be about twice 
larger than the boost factor in the toy model with 4 
electrons and 4 positrons in the final state. 

In Fig. [6l we compare the fiux from a subhalo with 
the fiux from the host halo 



1 



host 



Mb 



dEo k{pI^So)Q{Eo). (50) 



iTT b{E) Je 

Since the dependence on the host halo density pro- 
file is insignificant, we use the homogeneous DM dis- 
tribution with pdm = 0.3 GeVcm^^ and (crw)o — 
3 X lO^'^^cm'^s"^. We have chosen a local boost factor 
of 5*0 = 600 and a DM particle mass of Mdm = 2.1 TeV 
in order to fit the ATIC and PAMELA data (see Fig. 
[7]). Thus the fiux in the presence of significant contri- 
bution from a subhalo is consistent with the ATIC and 
PPB-BETS data but may contradict the Fermi data if 
we simultaneously fit to the PAMELA positron ratio. 

The flux from the host halo without subhalos is 
shown in Fig. [51 where we use the same model as above 
but increase the boost factor to Sq = 900 in order to fit 
the PAMELA, Fermi, and HESS points. For a smaller 
DM mass, one can also fit the PAMELA and ATIC data 
(as was shown, for instance, in [ist). Thus either the 
ATIC or the Fermi data can be explained by DM an- 
nihilation in the host halo for different parameters in 
the DM models (but not simultaneously, since ATIC 
and Fermi data have significant deviations from each 
other at energies between 300 GeV and 700 Gev) . The 
fits of the fiux from the host halo to the Fermi-LAT 
and PAMELA data and the corresponding ranges of 
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parameters for a more realistic counterpart of our toy 
model with the annihilation chain 2% 201 ^ 
can be found in [tH] • 

V. CONCLUSIONS 

In this paper we have studied the e+e^ flux from 
annihilating DM. We analyzed the dependence of the 
flux on the choice of DM model, on the profile of the 
DM halo, and on the presence of a local subhalo. 

The DM models we have considered here are charac- 
terized by the mass of the DM particle A/dm and by the 
number of step k in the annihilation-decay process from 
DM to e+e". In particular, we studied a toy model 
where the DM particles annihilate into two scalars that 
decay by a chain of (k - 1) two-body decays into elec- 
trons and positrons. In this model, log ~ 
and, as expected, F ~ for E <^ E^,. The behavior 
of the flux for E ^ E^, \s model dependent. The spec- 
trum has a sharper cutoff in models with fewer decay 
steps. 

The dependence on the shape of the DM host halo is 
very mild, provided that the DM distribution is suffi- 
ciently smooth. For the typical DM halo profiles, such 
as NFW, Einasto, and Isothermal, 5n < 10%. 

If a small number of DM subhalos contribute signifi- 
cantly to the DM annihilation, then the corresponding 
source function may have large variations and the in- 
dex of the propagated e^e^ flux will change as well. 
In general the index grows, i.e., the flux from a sub- 
halo is harder than the flux from the host halo. Using 
the Via Lactea II simulation, we argue that, on aver- 
age, we expect one M > IO^Mq subhalo within 3 kpc 
from the Earth, and the flux from such a subhalo will 
leave a significant imprint on the electron spectrum 
above ~ 100 GeV, if its distance is less than ^ 1 kpc. 
Thus, based on the Via Lactea //simulation, there is 
at least a ~ 4% chance to observe the flux of electrons 
and positrons from a local DM subhalo. Extrapolating 
below the VL2 resolution limit, we estimate that this 
probability may grow to ^15%. 

In the presence of DM subhalos, there may be some 
features in the e'^e~ flux spectrum at high energies, 
but at low energies we should still expect the univer- 



sal index n « —2 of the DM annihilation flux. Future 
observations will help to distinguish between the fol- 
lowing possibilities: 

• The additional flux is dominated by a local DM 
clump at all energies. Then the index of the flux 
is not universal and, as a rule, n > —2, i.e. the 
flux is harder than the flux from the host halo. 
However, as was pointed out in [15J, we need a 
flux with an index n = —2.2 ± 0.2 in order to 
fit both ATIC and PAMELA data. This prob- 
lem becomes even worse if one tries to fit simul- 
taneously the PAMELA points and rather soft 
Fermi/LAT spectrum. Therefore this possibility 
contradicts current data. 

• The additional flux is dominated by a local DM 
clump at high energies, 100 - 1000 GeV, and by 
the host halo at low energies, 10 - 100 GeV. We 
present an example with a DM clump at 1 kpc 
from the Earth consistent with the ATIC and 
PAMELA data. However, this flux may be in- 
consistent with the Fermi/LAT and the HESS 
data. The flux from annihilating DM has an in- 
dex n « —2 at energies E < 100 GeV and n > —2 
at energies 100 GeV < < 300 GeV. 

• The additional flux is dominated by the host halo 
at all energies. For different DM models this flux 
can fit either PAMELA and ATIC or PAMELA 
and Fermi data but not both. The property that 
the positron ratio is fitted in both cases is a con- 
sequence of a general fact that, independently of 
DM model and the shape of the DM halo profile, 
we expect an index n — 2 for energies much 
smaller than the cutoff scale E <^ E^,. 

Our main conclusion is that at lower energies E <C 
A/dm, i-e. in the PAMELA range 10 - 100 GeV for 
A/dm ^ 1 TeV, one should expect the flux from the 
DM to have a universal index tlk, —2. At higher ener- 
gies, i.e., in the ATIC - Fermi range 100 - 1000 GeV, the 
behavior of the flux is model dependent. In the pres- 
ence of a significant nearby clump of DM one should 
expect an ATIC-like bump in the spectrum whereas in 
the absence of large nearby clumps the flux should look 
smooth, similar to the Fermi/LAT and HESS data. 
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